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Abstract — This paper addresses the optimal design of limited- 
feedbacli downlinli multi-user spatial multiplexing systems. A 
multiple-antenna base-station is assumed to serve multiple single- 
antenna users, who quantize and feed back their channel state in- 
formation (CSI) through a shared rate-limited feedback channel. 
The optimization problem is cast in the form of minimizing the 
average transmission power at the base-station subject to users' 
target signal-to-interference-plus-noise ratios (SINR) and outage 
probability constraints. The goal is to derive the feedback bit 
allocations among the users and the corresponding channel mag- 
nitude and direction quantization codebooks in a high-resolution 
quantization regime. Toward this end, this paper develops an op- 
timization framework using approximate analytical closed-form 
solutions, the accuracy of which is then verified by numerical 
results. The results show that, for channels in the real space, the 
number of channel direction quantization bits should be (A/ 1) 
times the number of channel magnitude quantization bits, where 
M is the number of base-station antennas. Moreover, users with 
higher requested quality-of -service (QoS), i.e. lower target outage 
probabilities, and higher requested downlink rates, i.e. higher 
target SINR's, should use larger shares of the feedback rate. It 
is also shown that, for the target QoS parameters to be feasible, 
the total feedback bandwidth should scale logarithmically with 
the geometric mean of the target SINR values and the geometric 
mean of the inverse target outage probabilities. In particular, 
the minimum required feedback rate is shown to increase if 
the users' target parameters deviate from the corresponding 
geometric means. Finally, the paper shows that, as the total 
number of feedback bits B increases, the performance of the 
limited-feedback system approaches the perfect-CSI system as 

Index Terms — Beamforming, bit allocation, channel quanti- 
zation, limited feedback, multiple antennas, outage probability, 
power control, spatial multiplexing. 



I. Introduction 

Multiple-input multiple-output (MIMO) technology can po- 
tentially provide significant performance improvements for 
wireless systems. More specifically, in the context of cellular 
communications, the availability of multiple antennas at the 
base-station allows it to simultaneously transmit to multiple 
users by multiplexing their data streams and hence improve 
the total downlink rate. Such systems are generally referred to 
as multi-user spatial multiplexing systems. 

The performance of multi-user spatial multiplexing systems 
depends heavily on the amount of channel state information 
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(CSI) at the base-station |[T], Such information is neces- 
sary to form the downlink transmission beams at the base- 
station and to perform rate/power adaptation for each user. 
Acquiring channel information, however, is a challenging issue 
in practice especially for frequency-division duplex (FDD) 
systems, where the uplink and downlink channels use different 
frequency bands. In these systems, the users need to explicitly 
quantize their channels and send back the quantized informa- 
tion through a shared feedback channel. The feedback link in 
such a system is usually a rate-Umited control channel, hence 
the term limited-feedback systems. 

The scarcity of feedback bandwidth in limited-feedback 
systems necessitates efficient channel quantization codebook 
structures and an optimal allocation of feedback bits among 
those codebooks. This paper aims at deriving these structures 
and the corresponding bit allocation laws for limited-feedback 
multi-user spatial multiplexing systems. 

A. Related Work 

Multi-antenna communications with limited CSI is exten- 
sively studied in the literature for single-cell systems ITl- lfTOl 
and to some extent for cooperative multi-cell networks ifTTl - 
[13 J. This paper focuses on single-cell multi-user systems. 

The major advantage of multi-antenna multi-user systems 
with respect to single-user systems is the sum-rate multiplexing 
gain, which follows from the fact that multiple simultaneous 
transmissions can be established in a multi-antenna downlink. 
In order to preserve this gain in limited-feedback systems, the 
author of f2\ shows that the total feedback rate should scale 
linearly with the signal-to-noise ratio (SNR) in dB scale. The 
work in ||2l addresses a setup with small number of users. In a 
network with large number of users, there is another source of 
sum-rate improvement, referred to as multi-user diversity gain, 
which is realized by the base-station opportunistically schedul- 
ing users with favorable channel conditions. For scheduling, 
a well justified approach is to choose users with high channel 
gains and near-orthogonal channel directions H], I!]-!?), ifTOl . 
The authors of 13] specifically show that one needs the channel 
gain information (CGI) in addition to the quantized channel 
direction information (CDI) in order to realize the multi-user 
diversity gain. The gain information however is assumed to 
be perfect in pj. The split of feedback bits between CGI and 
CDI quantization introduces an interesting tradeoff between 
multiplexing gain and diversity gain. This tradeoff is studied 
by Q, where the authors numerically show that more bits 
should be used for CGI quantization in order to benefit from 
multi-user diversity gain as the number of users increases. 

For the purpose of precoding the information intended 
to scheduled users, two distinct approaches are proposed 
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in the literature. The first group of work uses zero-forcing 
beamforming based on the quantized directions [31, ||6l. The 
second group adopts a codebook (or muhiple codebooks) of 
orthonormal beamforming vectors and selects beamforming 
directions based on the signal-to-interference-plus-noise ratio 
(SINR) feedbacks from the users U, g), ITOI. The authors 
of IfTOl specifically claim that in the regime of large number 
of users, the orthogonal beamforming approach outperforms 
the zero-forcing method. Finally, the authors of jS] present a 
more thorough analysis of multi-user limited-feedback systems 
considering joint training, scheduling, and beamforming. 

B. System Model 

This paper considers a limited-feedback multi-user system 
with a Af -antenna base-station and M single-antenna users. 
The users send their quantized CSI to the base-station through 
a feedback link with a capacity of B bits per each downlink 
transmission block. Based on this quantized information, the 
base-station then comes up with the downlink transmission 
powers and beamforming vectors. 

Our goal is to optimize the system performance subject 
to the total feedback rate constraint. Although the proposed 
approach is rather generic, we impose several simplifying 
assumptions on the system model in order to achieve closed- 
form solutions. These assumptions are listed below and will be 
explicitly mentioned and justified whenever needed throughout 
the text: 
Assumptions: 

Al. Most of the analysis in this paper is based on the 
assumption of high resolution quantization, i.e. B — > oo. 
This assumption is made mainly to achieve a tractable 
formulation for system optimization and has been used 
frequently in the literature of quantization theory lfT4ll . 
ifTSll . In addition, the paper investigates, both numerically 
and analytically, the system performance with moderate 
values of B, and examines the regime where the high 
resolution results are applicable. 

A2. Users' channels are independent and identically dis- 
tributed (i.i.d.). 

A3. Users' channel directions are uniformly distributed over 
the A/-dimensional unit hypersphere. This includes the 
Rayleigh i.i.d. channels as an special case lfT6ll . 

A4. Users' channel magnitudes are independent from channel 
directions and can have an arbitrary distribution. 

A5. We assume a product structure for the channel quantiza- 
tion codebook, i.e. channel magnitudes and channel direc- 
tions are quantized independently. This product structure 
is justified in [Sl, ifTTl . According to the high resolution 
assumption, if we denote the magnitude and direction 
codebook sizes by Nk and Nk, we have Nk,Nk — > oo 
for each user 1 < fc < AL 

A6. The beamforming vectors are zero-forcing directions for 
the quantized directions. This assumption is to mimic 
the perfect-CSI case, where zero-forcing beamforming is 
shown to achieve asymptotically optimal sum-rate scaling 
with SNR 1T81. lfT9l. 

A7. The user channels are assumed to be real vectors. Al- 
though this assumption appears in the earlier literature. 



e.g. II20I . 1211 . we use it mainly for the ease of geometric 
representation of the quantization regions and the corre- 
sponding calculations. The extension of the analysis to 
complex space is discussed in Section [Villi 

C. Problem Formulation 

We formulate the system design problem as the minimiza- 
tion of the average sum power subject to the users' outage 
probability constraints. In order to differentiate between the 
users' QoS requirements, we assume different target SINR's 
and outage probabilities across the users. Our goal is to 
derive the optimal split of feedback bits among the users 
and the corresponding magnitude and direction quantization 
codebooks. 

Different variations of the power minimization formulation 
are used in the literature, e.g. in Il22l - ||26L as an appropriate 
formulation for fixed-rate delay-sensitive applications, e.g. 
voice over IP, video conferencing, and interactive gaming. 
The reliability of the fixed-rate link is usually achieved by 
applying power control at the base-station to compensate for 
the channel fading, as in Wideband Code Division Multiple 
Access (WCDMA) system standards ll27l . 

An alternative problem formulation is to maximize the 
average sum rate subject to a power constraint m-H, 0, 
llOl. In most of these formulations, the transmission power is 
fixed and the quantized information is used only to adapt the 
transmission rates. This type of formulation is more appropri- 
ate for variable-rate communication systems, e.g. Worldwide 
Interoperability for Microwave Access (WiMAX) and 3GPP 
Long Term Evolution (LTE) system standards ||28]| . 1291 . 

We make the following comments on the existing rate- 
maximization formulations in the literature: 

1. The sum-rate maximization problem assumes equal pri- 
ority for all users in terms of their applications. Since the 
users are not differentiated based on their QoS measures, 
this formulation cannot answer the question of how 
to optimally split the feedback bits among users with 
different QoS requirements. 

2. The channel gain information (CGI) is clearly an impor- 
tant factor in scheduling the users and also in setting the 
downlink rate for each user. However, most of the existing 
literature either assumes perfect channel gain information 
or completely ignores this information. It is therefore not 
clear how to optimally split the feedback bits between 
channel direction and channel gain quantizers in the 
context of sum-rate maximization problem. The works 
that look into this problem are mainly numerical and lack 
closed-form solutions Q, ITOl . 

One solution to the first issue raised above is to prioritize 
users with certain weights and consider the weighted sum-rate 
maximization instead of the sum-rate itself. These weights 
can be set by the scheduler based on users' QoS require- 
ments. The proportional fairness scheduler for example sets 
the weights based on users' backlogged traffic by assigning 
a higher weight to a user with larger backlog. This type of 
formulation appears for example in l30l - l33ll for perfect-CSI 
systems. Generalizing this formulation to limited-feedback 
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systems however appears to be a difficult problem. We are not 
specifically aware of any work that addresses such a problem. 
For this reason, this paper resorts to a power minimization 
formulation, which can easily incorporate the QoS constraints 
by assuming different target SINRs and outage probabilities 
across the users. The power minimization formulation also 
simplifies the CGI/CDI bit allocation problem and provides 
insight to the system design by allowing for closed-form 
asymptotic bit allocation solutions. 

D. Proposed Approach 

Our approach for solving the power minimization problem 
is to fix the channel outage regions in advance and transform 
the problem to a robust optimization problem. Assuming zero- 
forcing beamforming vectors, we first formulate the robust 
power control problem in the form of a semi-definite pro- 
gramming (SDP) problem. Using an approximate upper bound 
solution to the SDP problem, we then derive the codebook 
structures and the corresponding bit allocations in the asymp- 
totic regime where B — )• oo. Within the proposed approximate 
optimization framework, we show that for channels in real 
space: 

1. The optimal number of channel direction quantization 
bits is A/ — 1 times the number of channel magnitude 
quantization bits, where M is the number of base-station 
antennas. 

2. The share of the fcth user from the total feedback rate 
is controlled by log7fc and logl/qfc, where 7fe and 
qk are the user's target SINR and outage probability. 
As a general rule, a user with a lower target outage 
probability and higher target SINR needs a higher channel 
quantization resolution and therefore requires a larger 
share of the total feedback rate. 

3. For the outage probability constraints to be feasible, 
the total feedback rate should scale logarithmically with 
7, the geometric mean of the target SINR values, and 
\/q, the geometric mean of the inverse target outage 
probabilities. Moreover, the minimum required feedback 
rate increases if the users' target parameters deviate 
from the average parameters 7 and q, i.e. there is a 
feedback rate penalty for serving users with non-similar 
target parameters. The higher the deviation, the higher 
the penalty. 

4. As the total feedback rate B increases, the performance 
of the limited CSI system approaches the performance of 

_ B 

the perfect-CSI system as 2 TF^ . 
These optimality results are based on minimizing an upper 
bound of the sum power The closeness of the upper bound to 
the exact sum power is verified numerically in the paper 

E. Organization of the Paper 

The remainder of this paper is organized as follows. Section 
nil provides an overview of the perfect CSI system. Section 
nm presents the system design problem in its general form 
and describes our approach in transforming it to a robust 
design problem. Section |IV] describes the product channel 



quantization codebook structure. In Section |V] we study the 
power control optimization for fixed quantization codebooks 
and derive an upper bound for the average sum power. By 
using the sum power upper bound, we then optimize the 
product codebook structures in Section |VT] and derive the 
asymptotically optimal bit allocation laws. Finally, Section lVIll 
presents the numerical results and Section IVIIII concludes the 
paper. 

Notations: Most of the computations in this paper are in real 
space. The logarithm functions are base 2. The angle between 
any two unit-norm vectors u and v is defined as Z(u,v) = 
arccos|u^v| so that < Z(u,v) < 7r/2. 

II. Multiuser Spatial Multiplexing System with 
Perfect CSI: Outage is Inevitable 

We start by assuming perfect CSI at the base-station and 
show that, unlike a single-user system, outage is inevitable in 
the multi-user system even with perfect CSI. The difference 
with the single-user case is due to the fact that the base- 
station in a multi-user system needs to distinguish the users 
spatially and when the user channels are closely aligned it is 
not possible to satisfy the users' target SINR's with a bounded 
average transmission power. 

Consider a multi-user downlink channel with M antennas 
at the base-station and M users each with a single antenna. 
Let hfe e R^^, Vfe G R*^, Pk, and 7^ denote respectively, the 
user channel, the unit-norm beamforming vector, the allocated 
power, and the target SINR for the fcth user, 1 <k < M. The 
minimization of the transmission sum power subject the user 
SINR constraints is formulated as follows: 

M 



Pk 




|2 









>7fc, fc = l,2,---,M 



where the receiver noise power is assumed to be 1 for all users. 

A suboptimal solution for problem ([T]) is to use zero-forcing 
(ZF) beamforming vectors Vfc to eliminate the interference 
and find the power values Pk that satisfy the constraints with 
equality. This solution is asymptotically optimal in the high 
SNR regime ifTSl , |fT9l . Clearly, for the zero-forcing solution 
to be applicable, the users' vector channels need to be linearly 
independent. Since the channels are assumed to be independent 
random vectors, this condition is satisfied almost surely, i.e., 
with probability one. 

An important matter to consider with this solution is that the 
transmission powers would need to be extremely high when 
the users' channels are closely aligned, as the ZF beamforming 
vectors would be almost perpendicular to the corresponding 
channels in such cases. Therefore, it is not possible to always 
satisfy the SINR constraints with a bounded average power 
and as a result, a certain degree of outage must be tolerated 
by the users. 

To see this rigorously, define 6^ ~ Z(hfe,H_j.), where 
< fi'fc < f , and H_fc = span({h;|^ 7^ fc}). For zero-forcing 
beamforming vectors we have Z(vfc, hfe) = ^— 6'/j. Assume 
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that the users' channels are i.i.d. with uniformly distributed 
directions and independent channel magnitudes (with arbitrary 
distributions). The average sum power of the zero-forcing 
method is given by 



M 
k=l 



Ik 



M 
k=l 



7feE 



1 



E 



1 



As Ok is uniformly distributed in [0, the expectation of 
1/ sin'^(6'fc) becomes unbounded. 

To avoid unbounded transmit power, the users should tol- 
erate certain degrees of outage. A reasonable approach is to 
declare outage for user fc, i.e. set Pk = 0, when 

Q<ek< 01 

where 0^ ^ 1 is the smallest acceptable angle between and 
H_fe. With this assumption the average sum power is given 
by 



M 



fe=l 

M 
k=l 



1 

V/2 



tt/2 



dOt 



7feE[l/||hfcf] cot( 



2P, 



M 

Mu.csi TA; 



-1 "k 



k=l 



where the approximation holds for 6*^ <C 1 and 



def 



E[l/||hfc 



(2) 



(3) 



for i.i.d. users. The corresponding outage probabilities are 
Poui,k = 201/ TT for 1 < A: < M. 

Having studied the perfect CSI system, the next section 
describes a general framework for the limited-feedback system 
design. The insights achieved by studying this general form are 
used in the later sections for system design and optimization 
with product channel quantization codebooks. 

III. System Design Problem and Vector Channel 
Quantization: General Form 

To clarify the arguments, we start by some basic definitions. 
By a vector channel quantization codebook C of size N, we 
mean a partition of R*^ into N disjoint quantization regions 

S-^"), l<n<iV: 

For every quantization codebook C, we also define a quanti- 
zation function 

5(h) : R^^ ^ C, 

which returns the quantization region that h e belongs 
to. 

Now, for each user l<k<M, associate a codebook Ck of 
size Nk and the corresponding quantization function iSfe(hfc), 
where is the fcth user's channel. Further, define the ordered 
M -tuples 



M 



S(H) = [5i(hi),52(h2),--- ,5A,(hA/)] e n^*^- 



fe=i 



For a given total number of quantization (feedback) bits _B, 
target SINR values 7fe, and target outage probabilities g^, the 
system design problem is formulated as follows: 



mm 

Vfc(S(H)) 



H 



■ M 

E 

.k=l 



Pkim)) 



(4) 



M 



S.t. 



k=l 



prob 



P,.(S(H))|h^Vfc(S(H))| 
EP;(S(H))|h^v,(S(H))| 

l=jLk 



-<lk 



1,2, 



< 



qk, 



,M 



where the optimization is over the quantization codebooks Ck, 
codebook sizes Nk, the power control functions Pfc(S(H)) : 
n^ii Cfc ^ the beamforming functions Vfc(S(H)) : 

rife^i ilj\f, where is the unit hypersphere in R^^ . 

An exact solution to this problem is intractable. Our ap- 
proach in simplifying the problem is to fix the outage scenarios 
in advance and transform the design problem to a robust design 
problem that guarantees the target SINR's for the no-outage 
scenarios. 

Define the outage region ilk C Yik for user k such that 
prob[S(H) G flk] = Qk- Also define /fc(S(H)) as the activity 
flag for user k: 

4.(S(H))=I(S(H)el]^.), 

where !{■) is the logic true function. Whenever a user's 
channel resides outside the user's predefined outage region, 
the activity flag is on and the user must be served by the 
base-station, i.e. the user should not face an outage. 

Let us fix the codebook sizes Nk for now. For a robust 
system design, we need to design the codebooks, the power 
control functions, and the beamforming functions such that the 
target SINR's are guaranteed whenever /fe(S(H)) = l: 



mm 1 

Ck, 
Vfc(S(H)) 



M 



E^'^-(S(H)) 



.k=l 



S.t. 



inf 



Pfc(S(H))|w^Vfc(S(H))P 
we5.(h.) y: P/(S(H)) |w7^vKS(H))|Vi 

l^k 



(5) 



>7fe/fc(S(H)), 



V H e R'" and fc = 1,2,- - • ,Af 



(6) 



Note that by including the activity flag in the constraint (|6]l, 
this formulation guarantees the target SINR when /fe = 1 and 
returns Pk^O when lk~0- Also note that the activity flags are 
fixed in advance such that prob[//j(S(H))=0]=(7fc. 

The design problem in (|5]l is a complicated problem. In 
order to achieve a tractable reformulation, we accept two main 
simplifying assumptions (Assumptions A5 and A6 in Section 
EB: 

« We assume a product structure for the channel quan- 
tization codebook, where the channel magnitude and 
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the channel direction are quantized independently. Such 
a product structure, also known as shape-gain quanti- 
zation in the literature 1(141 . provides several practical 
advantages including faster quantization and lower stor- 
age requirement for the quantization codebooks. Product 
codebook structures are also shown to be sufficient struc- 
tures for effective interference management in multi-user 
systems fS). 

• We assume a fixed zero-forcing design for the following 
reason. Let us fix the quantization codebooks Ck and 
consider a time snapshot of the optimization problem in 
(|5]l with a fixed channel realization H and fixed quantized 
information S(H). For this specific point in time, we 
intend to optimize the beamforming vectors and power 
levels Pk such that the worst-case SINR conditions in 
^ are satisfied. It can be shown that this problem is 
non-convex with respect to Pk and and therefore 
difficult to solve in general. However, if one fixes the 
vectors v^, the problem would be convex in terms of 
Pk, since the objective and constraints are linear in Pk- 
Similar observations are made in 1341 - 11361 . In order to 
preserve the convexity structure, we therefore assume that 
the beamforming vectors Vfe(S(H)) are fixed as zero- 
forcing vectors for the quantized channel directions. This 
is to mimic the zero-forcing beamforming with perfect 
CSI, where the base-station knows the exact channel 
directions. The exact definition of the quantized channel 
directions is provided later in Section IIV-BI 
With these simplifying assumptions, the robust design prob- 
lem reduces to the following subproblems; 1) optimizing 
the power control function for fixed beamforming vectors 
and codebook structures; 2) optimizing the product codebook 
structure itself. The following sections address these subprob- 
lems. 

IV. Product Quantization Codebook Structure 

In this section, we describe the product quantization code- 
book structures and specify the corresponding outage regions. 
To be more exact, for a given target outage probability qk, we 
specify the magnitude and direction outage regions such that 

qk = qk + qk, 

where the magnitude outage probability qk is the probability 
that the channel magnitude resides in the specified magnitude 
outage region and the direction outage probability qk is the 
probability that the channel direction resides in the specified 
direction outage region. 

A. Magnitude Quantization Codebook and Magnitude Outage 
Region 

For each user 1 < fc < M, we use a magnitude quantization 
codebook 

for quantizing the channel magnitude squared Yk == Hh^lp. 

in) 

Here y^. are the quantization levels and Nk is the magnitude 
codebook size. 



For a given magnitude outage probability q, we define the 
magnitude outage region as the leftmost quantization interval 



[O, Uk j - The first quantization level is therefore fixed as 



F-\qk), 



(7) 



where F ^(•) is the inverse cumulative distribution function 

(cdf) ofYk = \\\ik\?.^ 

We further define Ck as the set of magnitude quantization 
regions, i.e. the set of quantization intervals for j|hfc|| = \/Yk: 



Ck 



t(1) t(2) 



where J^"'' 



J (Nk + l) def 

and yl ' = oo. 



(8) 



Note 



that the definition uses the square root of the levels as the 



quantization levels yj^ 
Finally, for Yk > 



(") 



are defined for quantizing ||hfc|| . 
we define the quantized magnitude 
Yk as the quantization level in that is in the immediate left 
of Yk, i.e.. 



Yk 



(ri) 

Vk 



if 



yl'<yk<yl ■ 



(9) 



For reasons that are clarified later in Section [VT] we are in- 
terested in a magnitude quantization codebook that minimizes 



E 



l/Yk 



It is shown in 1371 that the optimal codebook with 
such a criterion is uniform (in dB scale) in the asymptotic 
regime where Nk — > oo. We denote such an optimal codebook 
by Y*f. and refer to it as the uniform magnitude quantization 
codebook in the remainder of this paper. 

The work in l37l further shows that the uniform codebook 
satisfies the following upper bound: 



E 



1 

Y* 



< Pm 



1 + N, 



<k{Nk) 



(10) 



where 

p„„,„=E[i/y,] = E[i/|ih,f] 

as defined in On the left hand side of ( fTOl i. the variable 
Yj^ is the quantized magnitude variable associated with the 
uniform codebook Y^. On the right-hand side of ( fTOl ), 

nYk] 



def 
OJ = 



rfE[l/Yk] ' 



where 



7^= lim 

y— )-oo 



and /(•) is the probability density function (pdf) of Yk. The 
function (^k{n) depends on the magnitude outage probability 
(jk and is defined as the solution to the following equation: 

-1 



where y^^-* = F ^{(jk) as defined earlier. It can be shown that 



for any qk > and hence y^^-* > 0, we have 
lim Cfe(n) = 1- 



(11) 



The bound in ( flOb and the limit in ( fTTT ) are used in Sec- 
tion |Vl] for optimization of the product channel quantization 
codebook structure. 
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Fig. 1. Spherical cap B^"' around u^"' on the unit liyperspliere iXm- 



B. Direction Quantization Codebook and Direction Outage 
Region 

For each user l<fc<A/, we use a Grassmannian codebook 
Ufc of size Nk for direction quantization: 



{ui^),ur),..,uf^)}, 



(12) 



where u|."^ vectors are M-dimensional unit-norm Grassman- 
nian codewords. 

Every channel realization hj. is mapped to a vector 
UA:(hfc) e Ufe that has the smallest angle with h^: 



Ufe(hfc) = arg mill /(h^, u) 

ueUfc 



(13) 



The vector Ufe is referred to as the quantized direction for 
the channel realization hfc. The corresponding quantization 
regions, according to the Gilbert- Varshamov argument 
can be covered by the following spherical caps: 



Ck 



B 



(2) 
k ' ■ 



(14) 



where 



B 



") = {weUM|^(w,u("))<0,.} 



is the spherical cap around uj^"^ as shown in Fig. [T] In this 
definition, 

4)k = arcsin 4 

is the angular opening of the caps and 5k is the minimum 
chordal distance of Ufc. It should be noted that covering the 
direction quantization regions with the spherical caps enlarges 
the quantization regions. By considering the constraint (|6]l in 
the robust design problem (|5]l, such enlargement of the regions 
will lead to an upper bound for the average transmission 
power 

In order to describe the direction outage regions, define 



where 



U_fe = span({u(|Z ^ k}) 



and Ufc is the quantized direction for user k. This is a similar 
definition as in Section |ll] except that the exact channels hfe 
are replaced with the quantized directions Ufc. Similar to the 
discussion in Section for the perfect CSl system, we say 
that user k is in direction outage if 

Q<ek< ei, 

where 6'^ is the minimum acceptable angle between u^. and 
U_fc. This implicitly defines the direction outage regions of 



the users. Finally, assuming that Ok is approximately uniform 
in [0,7r/2j3, the direction outage probability is given by 



qk ~ -Wfe. 



(15) 



In Section IVII which addresses the product codebook 
optimization, we will need the following inequality, which 
describes the dependence between the angular opening (pk and 
the direction codebook size Nk'. 



sin (j)k < 4 A A/ TV, 



(16) 



where Xm^ {y^r{{AI+l) /2) /r(A//2)) "-^ . This inequality 
holds for large enough values of Nk and its proof is presented 
in ifJTl . The approximation on the left-hand side of (fT6b 
assumes (f> ^ 1, which is justified by the high resolution 
assumption Nk S> 1 (Assumption Al and A5 in Section IFbI) . 

C. Product Codebook Structure 

Using our definitions of the magnitude and direction quan- 
tization regions, we define product channel quantization code- 
book Ck for each user k as follows: 



Ck = {CkxCk)UOk, 



(17) 



where C and C are the magnitude and direction quantization 
regions in (O and (O, and 



Ok = <h 



is a ball centered at origin corresponding to the magnitude 
outage region. 

Based on the definitions of the magnitude and direction 
outage regions, the activity flag for user k is given by 



> 



(18) 



where I(-) is the logic true function, and 9k, 0'^, and y'lp are 
defined in Sections IIV-AI and IIV-BI 

According to the activity flag expression in ( fTSl ), the outage 
event for user k can be expressed as the union of the magnitude 
outage event, corresponding to the channel magnitude lying in 
the magnitude outage region, and the direction outage event, 
corresponding to the channel direction lying in the direction 
outage region. By using the union probability formula we have 

prob[/fe = 0] < 4fc + iik « F{y\^^) + -91 

TT 

where qk and qk are the magnitude and direction outage 
probabilities respectively and the approximation follows from 
( fTSl l. In order to satisfy the target outage probability qk, we 
therefore impose the following constrain^: 



F{v^^) + -9l<qk. 



(19) 



' This holds if the channel direction are uniformly distributed (Assumption 
A3 in Section irBt and the codebooks undergo sufficient random rotations. 

-Note that this constraint is stronger than prob[/fe = 0] < and would 
therefore lead to an upper bound on the objective function (average sum 
power). This agrees with the direction of our analysis in achieving an upper 
bound for the sum power 



7 




Fig. 2. Sector-type channel uncertainty region Sk = S{Rk, rf^,\ii^,4>k) for 
user k. 



Finally, the product channel quantization codebook size is 
given by 

\Ck\^Nk^NkNk + l, (20) 

where Nk and Nk are the magnitude and direction codebook 
sizes respectively. 

To summarize, for a given target outage probability g^. and 
a given codebook size N^., we propose a product quantization 
structure and specify its outage regions and the corresponding 
activity flags such that prob[/fe = 0] = qk- The proposed 
product codebook structure is parameterized by the magnitude 
and direction codebook sizes Nk and Nk, the minimum ac- 
ceptable channel magnitude and the minimum acceptable 
directional separation 6'^. 



V. Optimization of the Power Control Function 
WITH Sector-Type Channel Uncertainty Regions 

Assuming the product codebook structure in Section |IV] 
this section addresses the optimization of the power control 
function. For this purpose, we fix the quantization codebooks 
Cfe and the corresponding outage regions. Furthermore, as 
mentioned in Section |III] we make the simplifying assumption 
that the beamforming vectors are the zero-forcing vectors for 
the quantized directions. 

For the product quantization codebooks considered in this 
paper, the quantization (or channel uncertainty) regions are 
sector- type regions as shown in Fig. |2] A sector- type region 
is parameterized as 

S{R, r, u, 0) = I h e R'*^! < ||h|| < Vr, Z(h, u) < , 

where in the terminology of Section |IV] u is the quantized 
direction and r is the quantized magnitude, which is denoted 
as Y in 

For a specific point in time, consider the channel realizations 



T 

2 ' ■ 



H 

(or channel uncertainty) regions 



h^^] and the corresponding quantization 



Sk^S{Rk,rk,Uk,M=Sk{hk), 

where the quantization functions Sk{h.k), l<k<AI, are de- 
fined in Section |lll] Also, let Ik denote the corresponding 
activity flags. 

The goal is to optimize the power control function for the 
robust design problem (|5]l. Therefore, for the current channel 
realizations H, we have to find the transmission power levels 



Pk that minimize the instantaneous sum power subject to the 
worst-case SINR constraints: 

Al 

J2Pk (21) 



mm 

Pk 



k=l 



s.t. inf 



1 



>jklk, fc=l,2,---,M (22) 



l^k 

where the beamforming vectors are fixed, since the quan- 
tized directions Uj, are fixed. Let us refer to the users with 
/fc = as the silent users and the users with Ik = 1 as 
the active users and let the set /C denote the set of active 
users: K, = {l<k<M\Ik = 1}. In general /C is a random 
set depending on the channel realizations and the specified 
outage regions. Now, considering the power control problem 
in (l2Tl i. we note that if a user k is silent, i.e. Ik — 0, the 
corresponding SINR constraint in (|22] | is redundant as the 
problem returns Pk = for such a user. We therefore confine 
the SINR constraints in (|22] | to the set of active users, i.e. the 
indices k & IC. 

According to the robust design formulation in (|5]l, for any 
channel realization H = [h^ , , • • • , hj/], the base-station is 
required to serve (guarantee the target SINR's for) all active 
users k £ IC. The power control problem in (ISTT i therefore 
must be feasible for the active users. The following theorem 
presents a sufficient condition that guarantees feasibility. 

Theorem 1: To ensure the feasibility of the robust power 
control problem in (l2Tl i. it is sufficient to have the following 
for all 1< fc < M: 



M-l 



Nk > 4A 



sin arctan 



sm ( 



1 + ^{AI - l)7fe 



(23) 

Proof: See Appendix lAl ■ 
This condition is referred to as the minimum quantization 
codebook size (MQCS) condition in the remainder of this 
paper 

In the high resolution regime, where the codebook sizes 
tend to infinity, the MQCS conditions are clearly satisfied; 
therefore, feasibility of the power control problem is not an 
issue as far as the high resolution analysis is considered. This 
condition however plays a key role in finding the minimum 
number of feedback bits B for which the asymptotic bit 
allocation laws are applicable. This issue is discussed in further 
detail in Theorem |5] of Section IVll 

We are now ready to solve the power control problem in 
dSTt . An exact numerical solution to this problem can be 
obtained by transforming it into a semidefinite programming 
(SDP) problem as described in the following. 

Theorem 2: The problem in jTli is equivalent to the fol- 
lowing SDP problem for M > 3: 

min V Pk (24) 



s.t. 



— ^^fcVfeVfc PiViv'[h{Xk-Hk)'i-M- 



Hk 



Ik 



Xk > 



l^k 



COS," 



Mfc > 0, Pk > 0, fc G /C, 
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where Im is the MxM identity matrix. 

Proof: The proof is based on the Polyak's theorem in 
||39l . See Appendix |B] for details. ■ 

Although the SDP reformulation provides a numerically 
efficient solution, it does not give the minimized sum power in 
a closed form. The availability of a closed-form expression for 
sum power is crucial in optimizing the quantization codebook 
structures. Moreover, for practical systems where the base- 
station needs to continuously compute and update the trans- 
mission powers for each fading block, the SDP reformulation 
would be of a limited application. We therefore resort to a 
suboptimal solution to the problem in d^TT i that provides a 
closed-form upper bound to the sum power This upper bound 
solution is used later in Section |VT] as the objective function 
for optimization of the users' quantization codebooks. 

First, we bound the SINR terms as follows. For the sector- 
type regions Sk = S{Rk,rk,Uk,4'k), we have 

|2 



. , Pk Vfe 

mi ■ i 



(a) 



1 



>- 



l^k 

PkTk infweSfc Iw^vfep 



l^k 

(b) Pkr k sin"^ {0k- 



|W^V,|2- 



c) 



(25) 



(26) 



where w = w/||w||. The equality (a) holds since the SINR 
term is monotonic in ||w||, i.e. the minimum occurs on the 
spherical boundary region ||w||=y/rfc in Fig. |2] The equaUty 
(b) holds since Vfe's are the zero-forcing directions for Ufc's. 
To see this, consider 9k = ^(ufc, U-fc) as defined earlier By 
considering the zero-forcing principle, we have Z{vk,Uk) = 
^ — dk and Z(vi, Ufc) = for / 7^ k. Now, noting the definition 
of the sector- type region Sk, we have 

max Z(w, Vfe) = - - Ok + (t)k, 
7r 

mill Z(w, V;) = - - 
weSfr z 

By substituting the cosine of these angles in the numerator 
and denominator of dZSl l. we achieve the final expression in 



In order to obtain an upper bound on the sum power, we 
set the last term in ( l26l l to be equal to 7^: 

PkTk sin^ {Oi 



Ik- 



{Y^i^kPiW sin^ (t)k + l 

This is a set of linear equations in Pk, k E K,, where JC is the 
set of active users. By solving these equations and computing 
Pfc, we achieve the following upper bound for the sum 

power: 

'EkelC^k/l^k 



where ak 



P = V Pk - 

MU / y ^ k 

keic 

A sin" 
\^ 7fc sin" ipk 



1 - EfeGK "fc 



and /3k 



Tk sm 



(27) 



The 



subscript MU in P^^^ stands for multi-user. The closeness of 
the upper bound solution in (|27| i and the solution to the SDP 
problem in ( l24l i is verified numerically in Section IVIII 



The upper bound in jTTI i is a bound on the instantaneous 
sum power for a single snapshot of channel realizations in 
time. As the users' channels change over time, the quantized 
magnitudes = Yk, the quantized directions and the 
corresponding angles 9k = Z(ufe,U_fc) all change with 
time. The variables ak and /3k in (|27] | are therefore random 
variables. Since we are interested in the expected value of the 
sum power as the design objective in (|5]l, we use the following 
sum-power upper bound approximation so that the expectation 
operation can be applied conveniently. 

Theorem 3: In the asymptotic regime with large quantiza- 
tion codebook sizes and small values of (/>k, we have 



fce/c 



(28) 



kejc 



where 



ek^^il + Cl), fk 

Tk 



2jk 
Tk 



(Cfe+Cl), a=cot0fc. (29) 



Here, the notation g{(f>) = o(0), for an arbitrary function g{-), 
means that lim^^o g{4')/4' = 0- 

Proof: See Appendix |C] ■ 
So far, we have only considered the active users k G K,. 
In order to make the results applicable to the general case 
where some users might be in outage, we substitute jk with 
7fe//c in definitions of and fk in (l28T l so that users with 
Ik = contribute zero power to the sum-power upper bound. 
We therefore use the following expression for the average sum- 
power upper bound, where we have also replaced by the 
quantized magnitude Y^: 



e[Pmu]«EieN + E™'^'^- 



(30) 



where (k— cot 6k and and fk are redefined as follows: 



-(1+Cfcj, Jk - ~ {Qk + Cfcj- 



Yk 



Yk 



This concludes the optimization of the power control func- 
tion. In the next section, we use the average sum-power upper 
bound in ( [30] l to optimize the product quantization codebook 
structures and to derive the asymptotic bit allocation laws. 



VI. Product Codebook Optimization and 
Asymptotic Bit Allocation Laws 

In this section, we study the quantization codebook opti- 
mization. For this purpose, we use the average transmission 
power bound in ( [3Q] l in order to optimize the users' magni- 
tude and direction codebook sizes for a given feedback link 
capacity constraint and to derive the optimal bit allocation 
across the users and their magnitude and direction quantization 
codebooks. The optimization process is asymptotic in the 
feedback rate B and assumes large quantization codebook 
sizes, 7Vfe,iVfc > 1. 

Consider the sum-power upper bound in (|30] |. Assuming 
that 6k ~ Z(ufc,U_fe) is approximately uniform in [0,7r/2] 
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and using the definition of the activity flag in ( fTST l. we have 

r ~ 1 ^ r^"^ 2'Yz r 

E[efe] = 7feIE ^ / 1+Cfc <^0k « l/F,. 



E[/fc]=27fcE l/Y^ 



7r/2 
1 



7r/2 



27fc 



E 



where the approximations hold for ^'^^Cl. By substituting 
these in dSOl l. we achieve 



Af 



" 1 ■ 

















(31) 



The parameter 0^ in OTT i is controlled by the direction 
codebook size Nk as described by (fTSI l. The term E[l/Yfe] 
on the other hand is controlled by the magnitude codebook. 
The asymptotically optimal codebook that minimizes this term 
is the uniform (in dB) magnitude quantization codebook Y*. 
By setting Y = Y* and by using (fTOl ) and ( fT6b , we can bound 
the average sum power in (ISTT i as follows: 



M 

fe=i "fe 



-2CA,(Wfc) 

fe 



(32) 



where p„ucsi defined in (O and the variable u and the 
function Cfc(') depend on the magnitude outage probability 
and the distribution of the channel magnitude as described in 
Section HV^ 

Our goal is to minimize the average sum-power upper bound 
in ( [32] i in terms of the magnitude and direction quantization 
codebook parameters. The optimization constraints are as 
follows. Assuming a total number of feedback bits B, we have 
the following constraint on the codebook sizes: 



M 



(33) 



k=l 



The two other constraints are the target outage probability 
constraints given by il9[ and the MQCS conditions in ( |23] |. 
For the total feedback rate B, the target outage probabilities 
qk, and the target SINR values 7^, the product codebook 
optimization problem is therefore formulated as follows: 



M 

mm > — 

Nk,Nk rrf ^k 



M 



s.t. Jl [NkNk + 1 ) = 2^ 

k=l 

F{y^k^) + \ 



(34) 
(35) 



Nk> 4A7\/ / sin arctan 



l+v/(M-l)7fe 




In order to obtain a closed-form solution for the optimal 
product structure, we simplify this problem as follows. 

First, by assuming Nk,Nk ^ 1 and using the fact that 
limjY^^j^ Cfc(^fc) ~ 1, we use the following approximation 
for the objective function: 



1 + iV, 



-ujN 



1 + 



1 + N^' + 



1 + 



(38) 



We therefore have the following approximate upper bound: 



2Pm 



E 



4A 



M 



Nk 



(39) 



Next, we simplify the optimization constraints as follows. 
We approximate the first constraint (l35l l as Jlfcli ^k^k = 2^. 
Regarding the outage constraint in (|36] |, one can easily show 
that the objective function in ( |34] | is a decreasing function 
of and 9'^. The constraint in ( [36] l should therefore be 
satisfied with equality at the optimum. In order to simplify 
this constraint, we make the assumption that the magnitude 
and direction outage probabilities are equally likeljH : 



qk = F{yl ') = y 

2 qk 
Ik = -&k = y- 

(1) 



According to this assumption, y\ = F ^{qk/2) and 



(40) 
(41) 

(42) 



are fixed and the codebook optimization is only over the 
codebook sizes. Finally, since the optimization is asymptotic 
in the codebook sizes Nk and Nk, the last constraint in (|37] ) 
is redundant. This constraint however is used later to derive a 
lower bound on the total feedback rate B such that the target 
outage probabilities are feasible. 

Now, by using the approximation in ( [38] ), the optimization 
problem in (|34] | simplifies to the following optimization prob- 
lem: 



M 



Tfe 

M 

S.t. Y[NkNk = 2''. 



(43) 
(44) 



fe=i 



Define Bk ~ log Nk and Bk — log Nk as the number of 
magnitude and direction quantization bits respectively. 



^It can be shown that any other division of the form qf^ = aq^. and ijk = 
(1 — a)qf, with < a < 1 only changes the bit allocations in Theorem |4] 
by a finite constant and therefore does not affect the asymptotic bit allocation 
results. 
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■Bk (numerical) 
-Bk (analytical) 
■Bk (numerical 



k = 1 




100 



110 120 130 140 150 160 170 

Total number of quantization bits (B) 



180 



Fig. 3. Analytical magnitude and direction bit allocations in j47t and )48t 
vs. numerical bit allocations for three users with 71 = 15dB, qi = 0.02 and 
72 = 73 = lOdB and 92 = 93 = 0.05. 



Theorem 4: Define 



1 



B 



M2 



M-l, 1 
M ^ q 

M-l 1 
B + log ^ + K„ 



M 



(45) 
(46) 



where k., 



^log^riT^ and (rife is the ge- 

ometric mean of the target outage probabilities. The optimal 



values of B^ and are given by 



Bk — Ba 



lOf 



7fe 



Bfc =Save + (A/-l)l0g 



where 7= (Jlfc 7fc)^^'*^ is '^he geometric mean of the target 
SINR values. 

Proof: See Appendix ID] ■ 
Although the bit allocation laws in Theorem |4] are derived 
with the simplifying assumption of equal magnitude and 
direction outage probabilities in ( l40l i and dTTT i. the numerical 
results, as shown in Fig. |3] verify that the analytical results 
in (|47] | and (l48T l are close to the bit allocations derived by 
numerical minimization of ( l34l i. For the example in this figure, 
the base-station has A/ = 3 antennas and serves three users 
with the target parameters 7^. = 15dB and qk = 0.02 for 
the first user (k = 1) and jk = lOdB and q^ = 0.05 for 
the two other users (k = 2,3). The user channels are i.i.d. 
and hk ^ A/^(0, 1m), where Im is the MxM identity matrix. 
Also, the number of bits are rounded to the closest integer 
number^ 

Corollary 1: For each user k, the optimal number of mag- 
nitude and direction quantization bits are related as follows: 

Bk = {M-l)Bk + (M-l) log - + A/k„„, (49) 

Qk 



log — 

Qk 

Ik , 



7 



2(M-l)log- 



Qk 



(47) 
(48) 



This is a popular approach in solving integer programming problems, 
where the integer conditions are relaxed and the optimum of the relaxed 
problem is rounded to the closed integer 1401 . 



where is defined in Theorem|4] Moreover the total number 
of quantization bits for user k is given by 

Bk = Bk + Bk = ^B + A/log ^ + (2A/-1) log ^. (50) 
AI 7 qk 

As it is expected, if the users are homogenous in their 
requested target parameters, i.e. qk and jk are the same for 
all users, each user takes an equal share of jjB of the 
total feedback rate. In the case of heterogenous users, on the 
other hand, a user with a higher QoS (lower target outage 
probability) and a higher target downlink rate (higher target 
SINR) uses a higher feedback rate Bk- 

The bit allocation laws in Theorem |4] are asymptotic results 
in the feedback rate B ^ 00. In the following, in order to get 
a sense of how high the feedback rate should be, we determine 
a lower bound on B for which the target SINR values 7^ are 
in fact achievable with the target outage probabilities qk- 

Theorem 5: Assume 7^ > 1 and (7^ <C 1 and define 



Qk 



def -^Ik 

Qk 



For the target SINR's jk to be satisfied with outage probabil- 
ities qk, the following total feedback rate B is sufficient: 

B > iA/2log7+ (Af2-Af)logi + Af2logA + 6, (51) 
z q 

Here 7 and q are the geometric means of 7fe's and g^'s 
respectively and 

Q 



A 



min Qk ' 

i<fe<j\/ 



where Q is the geometric mean of Q^'s. The constant b in 
(ISTT l is defined as b^^AP+^AP\ogM+ApK^^j, where k„„ 
is defined in Theorem |4] 

Proof: See Appendix |E] ■ 

Several interesting results can be extracted from Theorem 
|5] First, we observe that for the QoS constraints to remain 
feasible, the system feedback link capacity should scale loga- 
rithmically with the geometric mean of the target SINR values 
and the geometric mean of the inverse target outage probabil- 
ities. Second, if we compare the case of homogenous users 
with the case of heterogenous users, we see that heterogenous 
users impose an additional requirement, AP log A, on the total 
feedback rate. If we think of Q^'s as users' QoS indicators, 
the variable A can be interpreted as a measure of discrepancy 
among users' QoS requirements. A higher QoS discrepancy 
requires a higher feedback bandwidth. 

Finally, in order to study the performance of the limited- 
feedback system as the feedback rate increases, we substitute 
the optimal magnitude and direction codebook sizes given by 
Theorem |4] into the average sum-power upper bound in ( [39l ). 
The following theorem shows the scaling of the average sum 
power with the feedback rate B. 

Theorem 6: For a limited-feedback system with a total 
number of B feedback bits, we have 



(52) 
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SDP solution 




60 65 70 75 80 85 90 95 
Direction quantization codebook size (N) 



100 



Fig. 4. SDP solution in j24j vs. 
users with target SINR's of 71 = 



the upper bound solution in )27t for three 
3dB and 72, 73 = 6dB. 
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Fig. 5. Bit allocations B/^ 
and 51=52=53=0.1. 



for three users with 71 = 



bits (B) 

2dB, 72=5dB, 73=8dB, 



where P^^csi is defined in (|2|i, q is the geometric mean of the 
target outage probabiUties, and 

_ 16M / 7r3/2(A/- l)r((A/ + l)/2) \ ^^^^ 
~ 7r(M-l) 16r(M/2) ) 

Proof: See Appendix |F] ■ 
If we define a quantization distortion measure as 

E fP ] - P 
D{B) = ^^^-^ (53) 

Theorem|6]implies that the distortion measure scales as 2^~m^ 
as £? — > 00. 

VII. Numerical Validation 

This section presents the numerical results that support and 
verify the analytical results in the earlier sections. 

1 ) Upper Bound Approximation for The SDP Problem in 
([24] |.- The codebook optimizations in Section [VT] are based on 
the sum-power upper bound solution in dZTj ) as an approxi- 
mation of the solution to the SDP problem in jTM . Here we 
investigate the accuracy of this approximation by comparing 
the two solutions for M = 3 users with target SINR's of 
7i = 3dB and 72, 73 = 6dB. 

To simplify the comparison, we assume perfect channel 
magnitude information, i.e. the quantized magnitude variables 
Tfc in (l24l i and jZ/j i are equal to the exact channel magnitudes. 
For direction quantization we use Grassmannian codebooks 
from BTI . The same codebook size is used for all users. The 
sum power values are averaged over 100 channel realizations 
for which the SDP problem is feasible. The user channels are 
i.i.d. and hfe --7V(0,Im)- 

Fig. |4] compares the two solutions as a function of the 
direction quantization codebook size N. As the figure shows, 
the two solutions converge as N increase^. This justifies the 

^This can mad e rigorous by showing that as TV^ — > 00 and (j)k ~^ 0> th^ 
inequahty in )25t is satisfied with an equality and therefore the upper bound 
solution is exact in the asymptotic high-resolution regime. 



use of the upper bound solution as the optimization objective 
in our high resolution analysis. 

2) Bit Allocations and Distortion Scaling: Fig. [3]in Section 
rvn compares the numerical and analytical bit allocations for 
M = 3 users with different target parameters. Here we repeat 
the process with a different set of parameters and record the 
share of each user Bk from the total number of feedback bits 
B. The target parameters are 7i=2dB, 72=5dB, 73=8dB, and 
(ji=(72=(Z3=0.1. Channel models are similar to those used in 
Figs. |3] and |4] and the bit allocations are rounded to the closest 
integer numbers. As the figure verifies, users with higher target 
SINR's receive larger shares of the total feedback rate. 

Finally, we investigate the system performance scaling with 
the number of feedback bits for the same set of parameters 
as in Fig. |5] For this purpose we use the average sum- 
power upper bound E [P^y] in ( |32] | and the definition of the 
distortion measure D{B) in i53[ . Fig. |6] shows the distortion 
measure as a function of B when numerical and analytical bit 
allocations are utilized. As expected, the two bit allocations 
show close performances. The figure also shows the distortion 
upper bound in Theorem |6] for the purpose of comparison. 

VIII. Concluding Remarks 

We conclude this paper by comparing the asymptotic 
magnitude-direction bit allocation law for the multi-user sys- 
tem with that of the single-user system discussed in 1371 . For 
the multi-user system and in the asymptotic regime where 
B — 00, the relation in (|49] | implies that the number of 
magnitude and direction quantization bits (for each user) are 
related as follows: 

4„ = (A/-1)4„. (54) 

The subscript MU in ( l54l i stands for multi-user. For single- 
user systems, on the other hand, we have the following bit 
allocation law |[37l : 

= — ^ — -B^jj, (55) 
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Fig. 6. Distortion measure with numerical and analytical bit allocations for 
three users with 7i=2dB, 72=5dB, 73=8dB, and <?i=g2='?3=0.1. 



with SU standing for single-user. If we define a relative 
quantization resolution as T — B/B, then we have 



MU ^ SU 1 



(56) 



which means that for the same number of magnitude quantiza- 
tion bits, the number of multi-user direction quantization bits 
is twice the number of single-user direction quantization bits. 
This is shown schematically in Fig. |7] As the figure implies, 
in order to make a single-user channel quantization codebook 
applicable to the multi-user system, each direction quantization 
region in the single-user codebook should be further quantized 
with the same resolution as the whole unit hypersphere. 

As the final note we mention that the results in this paper are 
based on real channel space assumption mainly for the ease of 
geometric representation of the quantization regions and the 
corresponding sum-power calculations in Section |V] However, 
the exact same approach introduced in this paper can be 
applied to complex space channels. The only main difference 
is in using the upper bound in ( fT6b for real Grassmannian 
codebooks. For complex Grassmannian codebooks we need to 
use the following bound [161: 



sin (hu < 2Ni. 



(57) 



By doing so, the rest of the analysis can be applied in a 
similar fashion to derive the bit allocation laws. In particular, 
the codebook optimization problem in (l43T l translates to the 
following problem for complex channels: 



M 

Elk 




k=l 'k 
M 



1 + n: 



2 •• 1 

— N 2(Af-l) 
k 



S.t. 



n ^kNk 



= 2' 



k=l 



By solving this problem, one can easily derive the bit alloca- 
tion laws and the system performance scaling for the complex 
space similar to the ones in Theorems |4] and |6] 

In particular, due to the difference between dTSI l and (|57] | 
in the exponents of the direction codebook sizes Nk, the 



asymptotic magnitude-direction bit allocation in the complex 
space turns out to be B = 2{M — 1)B for the multi-user case 
and B ~ {M — 1)B for the single-user case. Therefore, the 
complex-space quantization resolutions also satisfy the relation 
in ( l56b . 

Appendix 

A. Proof of Theorem Q] 

In order to prove the theorem, we use the following lemma: 
Lemma 1: To ensure the feasibility of the power control 
problem in ( |2TI ). it is sufficient to have 



tan 4>k 



< 



1 



(58) 



l + v/(A/-l)7fc' 

for all active users k G JC. 

Proof: The idea is to show that if the condition ( fSST i holds, 
the upper bound solution in ( l27b is a valid solution to the 
power control problem in ( 1211 1. For this purpose, it suffices to 
show that J2keK '^k<^, where is defined in (l27b . 
According to the condition ( fSSl ), we have 

1 



tan (pk 
sin 6k 



< 



1+ v/(A/-l)7fc' 



Then 

siii{9 k-(f>k) 
sin 0fe 

Therefore 



cos9k > - — ■ 1 > \/{M-l)-fk- 



tan (pk 



tancpk 



ak 



1 



sin^ {Ok - 



7fcsm (pk 

J2keK. ^k < 1, since |/C| < M. 



< 



M 



(59) 



(60) 



The MQCS condition in ( l23T l is equivalent to 



sin ( 



< sm arctan - , 

I I 1 + v/(M-l)7fc 



Combining this with the inequality in ( fT6b . we have 

• / / sine*" \\ 

sm(pk < sm arctan — , 

1^ \l + ^{M-lhkJJ 

which leads to 



tani/ifc < 



Slllf 



1 + v/(M-l)7fc' 



(61) 



(62) 



(63) 
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Now for all active users k Q K., we have 9k > 9^, since 6'jl 
is the smallest acceptable angle between user fc's channel and 
the subspace spanned by all other users' channels. Combining 
Ok > 01 with (I63]l leads to 



sm Ok 



1 + ^{M - l)7fc 

which according to Lemma [T] guarantees feasibility and there- 
fore completes the proof. 

B. Proof of Theorem |2] 

Consider the definition of the uncertainty region 

5fc = {w e R^^l < ||w|| </R, Z(w,u) <(j}Y 

As far as the constraint in (l22l i in concerned, the constraint 



||w|| < \/ R is redundant since the infimum occurs on the 
lower surface ||w|| = ^/r (see the process of deriving the 
upper bound solution in dZSt). 
Now define 



To 



Pk Iw^Vfcl 



1 



> Ik 



(64) 



E Pi |w^v, 
Ti - {w|||w|| > ^^}, 

T2 = {w ||w^Ufc| > ||w|| COS0} 

where the condition in ( |66] | is equivalent to Z(w, u) < 
With these definitions, the constraint in ( |22] | is equivalent to 



(65) 
(66) 



TinTz C To. 



(67) 



Now, define scalars t„ and symmetric matrices An, for n 
0,1,2, as follows: 



^0 = 
A2 

To = 



l^k 
I 



1 



COS^ (/ 
-1, Ti 



rUfcU^ 



0. 



With these definitions, the sets Tq. Ti, T2 can be expressed as 
sublevels of quadratic functions: 

Tn = {w |w'^yl„w < T„ } , n = 0, 1, 2. 

We therefore can use the following theorem from j39| to 
replace the condition ( |67] | with a SDP condition. This theorem 
is an important extension of what is known as S-procedure in 
the optimization literature |4()|. 

Theorem 7: Let M > 3 and An E R^'^^^'^ be symmetric 
matrices for = 0, 1, 2 and assume 

3 ui,i^2 e R s.t. i^iAi + V2A2 >~ 0. 

Define the quadratic functions /„(w) = w^y4„w. Then the 
following two statements are equivalent: 

I. /l(w) < Tl , /2(W) <T2 ^ /o(w) < To (68) 

Ao ^ XAi + fiA2 



II. 3 A > 0, /X > s.t. 



To > Ati + ^T2 



(69) 



For the problem in hand, one can easily find vi and 1^2 such 
that the condition uiAi +V2A2 is satisfied; therefore, the 
constraint in (|67] | translates to the SDP constraints in 

C. Proof of Theorem \3\ 

From the definitions of ak and jSk in dZTl i. we have 



7fcSin (pk 



7fc sin^ (t>k + sin^ [Ok - 4>k)' 



Ik/rk 



Pk 7fc sin^ (j)k + sin^ {Ok - 4>k) 
For small values of (pk <C 1, it is easy to verify that 

7fc sin^ (pk = •ykcpl + o(0fc), 

sin^ {9k - (f>k) = sin^ 61^ - {sm2ek)(pk + o{(j)k). 

After a few manipulations, one can show that 

ak = ekTkCpl +o{(j)l), 



Sk + fk4>k + 0{(j)k), 



and 



Pmv = -; = 2^ ek+fk<Pk+o[(j)k)- 



D. Proof of Theorem |4] 

By applying Lagrange multipliers method, we achieve 

1 Ik 



Nk = 



A 



(70) 



Nk 



1 



4Am 7fe 



AM-i ^M-l (6i°)2 



fl/-i 



(71) 



where A is a Lagrange multiplier that should satisfy the 
constraint HfcLi NkNk = N. By solving for A we get 

1 2J\/-1 



A = 



4A 



M 



M - 1 



n-j (ni 



(72) 

By substituting (|72l) in and (I2B, and using 9^ = jQk 
as in ( |42] |. and further manipulation, the optimal quantization 
resolutions = log Nk and = log Nk can be expressed 
as in (|47| and diST l. 

E. Proof of Theorem |5] 

For the target parameters to be feasible, the optimal direc- 
tion codebook sizes Nk—2-^'' are required to satisfy the MQCS 
conditions in i3l\ . With the assumption of qk <C 1, we have 
sin 6*^. w 6*^. < 1 and the MQCS conditions simplify to the 
following conditions: 



iv.>(^(i + y(I7^)) 



(73) 



In the following, we find a lower bound on B that guarantees 
the conditions in ( |73] |. 
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For AI > 2 and 7^ > 1, we have the following inequality: 

1 + ^W~Wk < v/2A/7fe. (74) 

To satisfy the conditions in ( |73] l, therefore, it is sufficient to 
satisfy 



A/-1 



By substituting TV^ from (fTTT i in ( fTSl l. we achieve 

X > (M-l)^V2M7fe. 
A 7fc 

Therefore it is sufficient to have 
A 



(75) 



(76) 



(77) 



By substituting the expression for A in ( |72] i into ( TTTl ). we 
achieve the following constraints on = 2^ for 1 < fc < M: 



TV > C- 7 



/7fe 



(78) 



where 7 = (0^7/.)'^'' and 9° = {U^kF" and 



/ 4A 



M{M~1) 



M - 1. 

Since jTSl ) is to be satisfied for all 1 < fc < A/, we have 
the following sufficient bound on = 2^: 



N>C- 7 



max — — 

l<k<M 



max - , ^ , 

l<fc<A/ 9° J 



(79) 



By substituting 91 = and ^° = f g in (l79]l and taking 
the logarithm of the both sides we achieve the lower bound in 
dSTl i. which completes the proof. 

F. Proof of Theorem |6] 

By substituting the optimal values of Nk and Nk given by 
(ITOI i and (ItTI i into the average sum-power upper bound in i 
we have 



2a 



M 

MU,CSI ^ Ik 



2Pm 



A/ 

E 



7fc 




+ A + (A/ - 1)A 



A/^A 



(a) -^Pmu.csi 



+ ^^'X7 ( ^ 



1 



where 7= (11^ Ikf'"' and 6*°= (Hfe 61°)''"' . For the equahty 
(a) we have used the expression (|72| i for A with N = 2^ and 

X=(4Am/(A/-1))(^'^-^'/^^ 



(80) 



By further manipulating dSOl l. we have 



2 ^ 




(81) 



where in (b) we use the fact the geometric mean is smaller 
than arithmetic mean: 



M 



M 

1 Ik 



and in deriving (c), we use the definition of P^ycs, in (|2|i and 
the fact that 9° ~ jq < 1, since ^ is a probability measure. 
By substituting 9° = f 9 in dSTT l. we achieve the bound in 
and the proof is complete. 
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